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TRIGONOMETRIC PARALLAX OF W51 MAIN/SOUTH 

M. Sato 1 ' 2 - 3 , M. J. Reid 2 , A. Brunthaler 4 , and K. M. Menten 4 

ABSTRACT 

We report measurement of the trigonometric parallax of W51 Main/South 
using the Very Long Baseline Array (VLBA). We measure a value of 
0.185±0.010 mas, corresponding to a distance of 5.4ll.o*|g kpc. W51 Main/South 
is a well-known massive star-forming region near the tangent point of the Sagit- 
tarius spiral arm of the Milky Way. Our distance to W51 yields an estimate of 
the distance to the Galactic center of Ro = 8.3 ± 0.46 (statistical) ±1.0 (sys- 
tematic) kpc by simple geometry. Combining the parallax and proper motion 
measurements for W51, we obtained the full-space motion of this massive star 
forming region. We find W51 is in a nearly circular orbit about the Galactic 
center. The H2O masers used for our parallax measurements trace four powerful 
bipolar outflows within a 0.4 pc size region, some of which are associated with 
dusty molecular hot cores and/or hyper- or ultra-compact H II regions. 

Subject headings: astrometry - Galaxy: fundamental parameters - Galaxy: kine- 
matics and dynamics - Galaxy: structure - ISM: individual (W51) - masers 

1. Introduction 

We are using the National Radio Astronomy Observatory'^ Very Long Baseline Array 
(VLBA), the Japanese VERA (VLBI Exploration of Radio Astrometry) array and the Euro- 
pean VLBI Network (EVN) to measure trigonometric parallaxes of maser sources in massive 
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The primary goal of these efforts is to delineate the spiral 



structure of the Milky Way by direct and accurate distance determination. Our current view 
of Galactic structure relies mostly on uncertain kinematic dist ances, which tend to overesti 



mate the source distances, sometimes over a factor of 2 (e.g., iReid et al.l l2009ct ISato et al. 



20101 ). Therefore, a precise map of the Galaxy based on trigonometric parallaxes of massive 
star-forming regions will significantly refine our knowledge of Galactic structure. In addition 
to parallax measurements, we also obtain source proper motions. Combining the proper 
motion and distance of the source with its radial velocity yields the full three-dimensional 
(3D) space motion of the source, allowing detailed study of the kinematics of the Galaxy. 

In this paper, we report parallax measurements of W51 Main/South, a well-studied 
massive star-forming region that hosts some of the strongest H2O mase r sources in the 
Gala xy. W51 Main/South is associated with a giant H II region W51 A ( 
2008b , where ongoing f ormation of massive stars has been reported (e.g., 



2000; 



Clark et al. 



'igueredo et al. 



Okumura et al. 



20091 ) . A reliable distance to W51 is important for the physics of massive 
star formation, since masses, luminosities and ages of young stellar objects depend sensitively 
on distance. Also, W51 has special importance for the study of Galactic structure, owing to 
its location very close to the tangency of the Sagittarius spiral arm, hence yielding a good 
estimate of the distance to the Galactic center by simple geometry. 



2. Observations and Data Analysis 
2.1. Observations 

We used the NRAO VLB A to observe the 22 GHz H2O maser source in the massive 
star-forming region W51 Main/South. The observations were conducted under VLBA pro- 
gram BR134, which included observations of four star-forming regions thought to be in the 
Sagittarius spiral arm of the Milky Way. Observations of the other sources will be reported 
in a forthcoming paper. 

Prior to the VLBA parallax observations, we used the NRAO Very Large Array (VLA; 
program AR677) in the most extended A configuration on 2008 September 28 to find back- 
ground extragalactic sources as p osition references n ear the target H 2 maser sources. Fol- 
lowing the methods described by IReid et al.l (j2009al ). we selected unresolved sources in the 
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NRAO VLA Sky Survey (NVSS) catalog within « 2° of each maser target and observed both 
background candidates and target masers at 8 and 22 GHz. For W51 Main/South, we found 
one new background source J1922+1504 for VLBA obser vations in addition to two known 



calibrators from the VLBA Calibrator Survey-1 (VCS1; iBeaslev et all 120021 ). J1922+1530 
and J1924+1540. However, at 0.3 mas resolution of the VLBA, we found the new source 
J1922+1504 displayed a resolved structure and was not useful for parallax measurements (see 
§3.ip . We also obtained an accurate position of the H 2 maser source in W51 Main/South 
to observe with the VLBA. 

The VLBA observations of W51 Main/South were scheduled on 2008 October 22, 2009 
April 27 and 30, and 2009 October 20, as listed in Table [H These dates well sample the peaks 
of the sinusoidal parallax signature in right ascension, maximizing the sensitivity of parallax 
detection and ensuring that the parallax and proper motion are essentially uncorrelated. We 
chose four optimal epochs to sample the parallax peaks in right ascension only, instead of 
sampling the peaks both in right ascension and declination over six or more epochs, since 
the parallax signature of W51 in right ascension is about twice as large as in declination. 
The two middle epochs were chosen close to each other and near one parallax peak (which 
we call a "1-2-1" schedule). The two middle epochs are needed for symmetry and to obtain 
similar accuracy for the parallax maximum and minimum. This yields zero correlation 
among the parallax and proper motion parameters. Compared to a "1-1-1-1" schedule over 
1.5 cycles (years) of the parallax signature, this "1-2-1" schedule also reduces the time span 
of observations to one year, reducing sensitivity to maser spot lifetimes. 



Our 22 GHz observations were set up in the same manner as described by iReid et al. 



(j2009a[ ) for 12 GHz CH3OH maser observations. For atmospheric calibrati on, four "geode 



tic" blocks of I CRF (the International Celestial Reference Frame) quasars (IMa et al. 



1998 



Fey et al.l l2004f ) with positions accurate to < 1 mas were observed at different antenna el- 
evations across the sky in left-circular polarization with eight 4 MHz bands that spanned 
450 MHz of bandwidth between 22.02 and 22.47 GHz and uniformly sampled all frequency 
differences. Using a four-block setup allows the middle block of observing time to be avail- 
able for phase-referencing rapid-switching scans between the target maser and background 
continuum sources when the source elevation was the highest at most stations. We measured 
multi-band delays and fringe rates, mostly due to un-modeled atmospheric propagation de- 
lays, to determine zenith delay errors as a function of time for each antenna. 

Table [2] lists source information for the target H 2 maser emission in W51 Main/South 
and background continuum sources. We placed three rapid-switching (phase-referencing) 
blocks between the geodetic blocks and switched every 30 s between the target maser and 
one of the three background continuum sources. For rapid-switching blocks, we employed four 
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adjacent frequency bands of 8 MHz bandwidth in both right- and left-circular polarizations. 
The four bands were centered at local standard of rest (LSR) velocities, Vlsr, of 276.0, 
168.0, 6 0.0, and —48 . km s " 1 , covering the maser velocity extent reported in previous 
studies (IGenzel et al.l Il981bl ; llmai et al.l 120021 ); H 2 maser emission was detected in the 
second, third, and fourth bands. 

The data correlation was performed at the VLBA correlation facility in Socorro, NM. 
The data from each antenna pair were cross-correlated with an integration time of 0.92 s. 
The expected time-averaging effect in wide-field mapping from this integration time, for a 
beam size of ~ 0.3 mas, is 5% lo ss in peak amplitude a t an angular distance of A0 « 2" 



from the correlation center (e.g., iThompson et al.ll200ll ). Each of the frequency bands was 



split into 256 spectral channels, yielding frequency and velocity r esolutions of 31.25 k Hz and 
0.42 km s -1 , respectively, for a rest frequency of 22.235080 GHz fjPickett et al.lll998l ) for the 
H 2 6i6 — >§23 transition. 

Figure [1] shows the 22 GHz H 2 maser spectrum of W51 Main/ South over all four 8 
MHz frequency bands from the first-epoch's data produced by scalar averaging of the data 
over the whole observing time at all antennas except Hancock (flagged) and Saint Croix (no 
data). 



2.2. Calibration 



We calibra ted the correla ted data using the NRAO Astronomic al Image Processin g 
System (AIPS; iGreisen 1120031 ) in the manner described in four steps by lReid et al.l (j2009al ). 
We performed "manual phase calibration" using J2253+1608 data to remove instrumental 
delay and phase and delay offsets among the four frequency bands. We then used one strong 
maser channel as the phase reference and subtracted reference phases of this channel from 
all data, including all maser channels and all background sources, by interpolating phase 
solutions in time. Since the position of J1924+1 540 was known to 0.3 mas accuracy (and 
that of J1922+1530 to 0.8 mas accuracy) (IBeasley et al.l 120021 ). the apparent map offset 
of J1924+1540 when phase referenced to a maser spot is due to a positional offset of the 
reference maser spot from the correlation center. We corrected for this positional offset by 
shifting the phase center of the reference maser spot with the AIPS task CLCOR and re- 
fitted the phase reference data. For the first through third epochs, we used the maser spot at 
Vlsr = 62.9 km s -1 in W51 South (feature 32 in Table[3]) as phase reference, after shifting the 

phase center by ( Act cos 5, AS) = (0''8816, K'1012) to the reference spot position. Since 

this spot became significantly weaker (from 754 to 21 Jy beam -1 ) and stronger emission 
appeared elsewhere in the same spectral channel as our observations progressed, we could 
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no longer use this spot as a phase reference at the fourth epoch and used another spot at 
Vlsr — 41.5 km s _1 in W51 South (feature 26 in Table [3]), after shifting the phase center by 
(Aacos<5, AS) = (K'2451, - - 1"1803) to its position. 

We flagged jumps greater than 60° in the reference phase in adjacent maser scans and 
all data between these scans. We also flagged all data from antennas that were not of good 
enough quality to give opacity correction and produced discrepant amplitude results. Tabled] 
lists all stations that were either not used during the observation or flagged in the calibration. 



2.3. Finding Maser Spots 

Since maser spots are spread over a large angle on the sky, we first made large-field 
image cubes involving all spectral channels for each 8 MHz frequency band from the first- 
epoch data covering a region of ~ 20" x 20" around the correlation phase center (i.e., the 
phase center position was shifted back with CLCOR from the reference maser position). For 
the large-field imaging, we used only shorter baselines among the inner five VLBA antennas: 
Fort Davis, Kitt Peak, Los Alamos, Owens Valley, and Pie T own. Robust weighting was 



adopted for weighting of the visibility data ( IBriggs et al.lll999l ). Each image cube was 4096 



pixels by 4096 pixels over 256 channels with pixel size of 1 mas, yielding a field of view 
of ~ 4" x 4". We mapped 27 sub-fields per channel: a 5 x 5 grid spaced by 4"096 plus 2 
sub-fields offset north and south of the center by 12"288. In order to remove the effects of 
sidelobes of strong maser features in the same spectral channel, all 27 sub-fields per channel 
were simultaneously CLEANed using the task IMAGR. 

After the 20" x 20" large-field imaging, we used Gaussian fitting with the AIPS task SAD 
(search and destroy) to detect maser spots. Our criteria for maser spot detection required 
that each detected spot had a peak intensity greater than seven times the rms noise level (7a) 
over two or more spectral channels within ±2 channels (i.e., within AVlsr =0.84 km s _1 ) 
at the same position within ±1 mas in each coordinate. From the first-epoch's data, we 
detected 1362 spots in 280 unique "features" (i.e., spots at the same position but in adjacent 
velocity channels). 



2.4. Mapping Maser Spot Distributions 

After the maser spot detection from large-field imaging, we imaged a small field of 1024 
pixels by 1024 pixels with pixel size 0.03 mas (~ 1/10 of the FWHM of the beam) centered 
at each spot position over detected channels, which yielded a field of view of ~ 30 mas x 
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30 mas. If two or more spots were in overlapping fields, we used one field of 2048 pixels by 
2048 pixels centered at the average position of the spots. All fields in the same channel were 
simultaneously CLEANed to remove sidelobes. We shifted phase centers using CLCOR to 
each field before mapping the field. (For a technical detail, since we could only shift the 
phase center to one position at a time, we shifted the phase center using CLCOR to one 
of the fields that were simultaneously CLEANed, then kept only the image of the phase- 
centered field by discarding the other field images, and repeated shifting and imaging for all 
fields.) We then mapped each field at all four epochs in exactly the same manner, except 
we doubled the field size for the fourth epoch to allow for a large proper motion of the spot 
and/or small residual position offsets of the field centers due to the change of the reference 
maser spot at the fourth epoch. 

After obtaining all four-epoch maps around each maser spot, we measured positions 
with the AIPS task SAD. We calculated spot positions relative to the reference maser spot 
at Vlsr = 62.9 km s -1 in W51 South (feature 32 in Table [3]), by subtracting the position 
of this spot from all spot positions. We then fitted a linear (relative) proper motion for 
each spot using the position versus time measurements. We identified a maser spot to be 
the "same" spot at different epochs if (1) the spot was detected at all four epochs as the 
strongest spot in the same imaging field, with an intensity greater than 0.5 Jy beam -1 and 
also greater than 10 times the rms noise of the map; (2) the fitted relative proper motion was 
smaller than 10 mas yr -1 (« 250 km s -1 at a distance of 5.4 kpc); and (3) the uncertainty 
for the fitted relative proper motion was smaller than 0.5 mas yr -1 . 

We grouped the identified maser spots into a maser feature (i.e., a group of spots with 
similar position and velocity) if (1) the positions were within 1 mas in each coordinate; (2) the 
radial velocities, Vlsr, agreed wit hin 7 km s" 1 ( which allows for maser emission from different 



hyperfine split components; e.g., [Walker 1119841 ); and (3) the fitted relative proper motions 



agreed within 0.3 mas yr -1 . Table [3] lists all 37 features with proper motions detected from 
all four epochs. We estimated the radial velocity for each feature as the intensity-weighted 
mean of the radial velocities of the spots. The position and proper motion of each feature 
relative to the reference spot were estimated as error-weighted mean of the values for the 
included spots. Figure [2] shows the proper motions of the detected features relative to the 
calculated mean motion of all regions (see Sections 13.21 and | 
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3. Results 
3.1. Parallax 

Since relative positions for two sources (i.e., maser and quasa r sources) at low antenna 



eleva tions are very sensitive to atmospheric delay errors (e.g., iHonma. Tamura fc Reid 



20081 ). we adopted an elevation cutoff of 30° for parallax measurement, below which data were 
discarded for each antenna. At each epoch, we imaged the background continuum sources 
after calibration (Figure [3]) and then fitted elliptical Gaussian brightness distributions to 
the sources using the task JMFIT. These positions were subtracted from the positions of 
the maser spot at Vlsr = 62.9 km s _1 in W51 South (feature 32 in Table [3]). Since this 
maser spot was the phase reference at the first through third epochs, it had essentially zero 
position. However, at the fourth epoch, this spot was not at the map center since the phase 
reference was a different maser spot. Since both the maser spot at Vlsr — 62.9 km s _1 and 
the continuum sources are equally affected by the change in reference feature, differencing 
their measured positions removes the effects of a different phase reference spot in the last 
epoch. 

The change in position of the maser spot at Vlsr = 62.9 km s -1 relative to background 
continuum sources was then modeled by the parallax sinusoid in east-west (Aa cos S) and 
north-south (5) directions and a linear proper motion. This requires five parameters, i.e., a 
single parallax parameter and two parameters for the proper motion in each coordinate. 
Independent "error floors" were added to positions in each coordinate to allow for unknown 
system atic errors and the n adjusted iteratively to make the reduced x 2 ^ 1.0 in each coor- 



dinate flReid et al.ll2009al ). 



Figure [3] shows contour maps of the background sources. Since one of the three back- 
ground continuum sources, J1922+1504, displayed heavily resolved structure (~ 0.3 mas), 
we did not use this source for parallax measurement. For the two quasars, J1922+1530 and 
J 1924+ 1540, after individual parallax fitting, we also attempted a "combined" fitting with 
the data from both quasars, solving for seven parameters, i.e., a common parallax and proper 
motion of the maser source against both background quasars but different position offsets 
for each quasar. 

Figure H] shows position measurements and parallax fits of the masers relative to the 
two background sources, J1922+1530 and J1924+1540. The results of parallax fitting are 
listed in Table |U We obtained the parallax of W51 Main/South to be it = 0.183 ±0.006 mas 
using J1922+1530 and vr = 0.187 ± 0.009 mas using J1924+1540. From the combined fit, 
we obtained 7r = 0.185 ± 0.007 mas. The uncertainty of each parallax fit was obtained 
from the formal fitting uncertainty. The adjusted error floors which yielded x 2 ~ 1-0 were 
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a x = 0.014 mas in the east-west direction and a y = 0.023 mas in the north-south direction 
using J1922+1530, and a x = 0.020 mas and a y = 0.037 mas using J1924+1540. For the 
combined fit, a x = 0.020 mas and a y = 0.038 mas. 

We also attempted parallax fitting using other maser spots in both W51 Main and 
South, which yielded essentially the same parallax as the maser spot at Vlsr = 62.9 km s -1 
in W51 South. This indicates that structural changes in the maser spots are unlikely to be 
a significant error source for our parallax measurement of W51 Main/South. The dominant 
error source of the position measurements is likely to be residuals in the atmospheric delay 
modeling, but structural changes in background quasars may also add some "noise." Random 
errors due to thermal noise are negligibly small, considering the high signal-to-noise ratio in 
both the maser and quasar maps. 

The error-weighted mean of the two measurements using each quasar is 0.184±0.005 mas. 
The deviation of each quasar measurement from the mean, < 0.003 mas gives an estimate of 
the error intrinsic to each quasar, including potential quasar structural changes and atmo- 
spheric delay error differences between the maser and quasar positions. In the right panel 
of Figure HJ the positions in the east-west direction at the second and third epochs spaced 
by 3 days differ by w 0.02 mas from each other, while the measurements on the same day 
using different quasars differ by a much smaller amount « 0.003 mas. The position difference 
between 3 days is likely due to the atmospheric delay errors, since structural change of the 
quasars over 3 days is unlikely. Therefore, the dominant source of the position errors appears 
to be the residuals in the atmospheric delay modeling. 

Since the two quasars lie relatively near each other on the sky (see Table [2] for angular 
separations and position angles relative to W51 Main/South), the atmospheric delay errors 
may be partially correlated for the two quasars. The uncertainty of the combined-fit parallax, 
estimated from the formal fitting uncertainty, assumes that the errors of the measurements 
using two quasars are uncorrelated. If they are 100% correlated, the estimated error needs 
to be multiplied by a/ (16 — 7)/(8 — 5) = a/3 (i.e., 16 data points and 7 parameters yielding 
9 degrees of freedom for uncorrelated data, and 8 data points and 5 parameters for a single 
quasar case). As an estimate of partially correlated errors, we average the factors of 1 
(for no correlation) and v3 (for a 100% correlation), and multiply the formal error by 
(1 + v / 3)/2 ~ 1.37, which yields an uncertainty of a n = 0.010 mas, as the best estimate for 
partially correlated data. We thus adopt the combined-fit parallax of n = 0.185 ±0.010 mas 
as the best estimate, which corresponds to a source distance of 5.4l1l '|g kpc. 

Our par a llax d istance is in good agreement with 5.1^4 kpc for W51 IRS 2 (or North) 



by Ku et al.l (120091) from 12 G Hz CH3OH maser parallax, but has much higher accuracy. 



The parallax by lXu et al.l ( 120091 ) has much larger uncertainty because the 12 GHz masers in 
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W51 IRS 2 had large spot sizes and were resolved out at longer baselines of the VLB A, which 
resulted in lack of data at the highest resolution. Also, the low signal-to-noise ratio due to the 
low 12 GHz maser flux density may have contributed to the error of the measurement. Previ- 



ous to the parallax measurements, the distance to W51 was mostly bas ed on statistic a 
laxes based on 22 GHz H 2 maser p roper motions, e .g., 7 ±1.5 kpc by iGenzel et all ( 



paral - 
1981bh 



for W 51 Main and 8.3±2.5 kpc by ISchneps et al.l (119811 ) and 6.1±1.3 kpc by llmai et al. 
fcpQ2h for W51 North. W51 Main/South and W51 IRS 2/North are separated by « l'.O east 
and 0'6 south. Both are thought to be associated with the H II region W 51 A. Our distance 
is consistent with the upper limit of 5.8 kpc by iBarbosa et al.l ( 120081 ) from spectroscopic 
study of W51 IRS 2, bu t is in significant disagr eement with the smaller distance of 2.0±0.3 
kpc recently reported by IFigueredo et al.l (120081) from spectrosc opic parallaxes of four O-type 
stars near the giant H II region W51 A. IFigueredo et al.l ( 120081 ) discuss an extra uncertainty 
in the derived distance that comes from the interstellar extinction law, however a factor of 
1.28 difference in distance they cite for adopting differe nt methods t o estim ate Ak will not 
account for the distance discrep ancy. As discussed by IClark et al.l ( 120091 ). the underesti- 
mated spectroscopic distance by IFigueredo et al.l ( 120081 ) may come from an underestimate 
of the luminosity class of the stars (04 to 07.5) used to derive the distance. Our distance 
of 5.4 kpc sug gests that the lum inosity of the P Cyg supergiant type star [OMN2000] LSI, 
whose nature I Clark et al.l ( 120091 ) discuss for three different assumptions on its distance, is 
indeed very high and of order ~5 x 10 5 L & . 



3.2. Maser Map and Internal Motions 



We detected 37 H 2 maser features in W51 Main/South that we could identify at all 
four epochs, and we list their relative proper motions in Table [3J We classified features into 
four separate groups according to their likely excitin g sources: hyper-compact (HC) or ultra- 
compact (UC) H II regions W51e2-NW, W51e2-E Jshi. Zhao fc Han IboiOaB ) and W51e8 
( IZhang fc Ho II1997I ). and a region we call "u" without detected millimeter or centimeter 
continuum emission (see Figured]). Note that a proper motion of 2 mas yr _1 corresponds to 
51 km s _1 at the distance of 5.4 kpc. 

Figure [2] gives a map of positions and internal motions of the detected maser features. 
The origin of the map is the position of the reference feature 32. The plotted internal motions 
are relative to the average velocity of the four regions, (Av x ) = —3 km s _1 and (Av y ) = 1 
km s _1 relative to the maser feature 32 (see Section I4.2p . The positions of the HC and 
UC H II regions are also plotted by crosses whose sizes approximate the beam sizes from 
previous observations: W51el, e3 and e4 by iGaume. Johnston fc Wilson I (119931 ): W51e8 by 
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Zhang fe Ho I (119971 ) and W51e2-W, e2-E, e2-NW and e2-N bv IShi. Zhao k Han I fl2010al ). 
Molecular-line studies sho w evidence for gravitational collapse of W51e2 (W51 Main) and 
W51e8 (W51 South) cores JHo fc Young Ill996l : lzhang. Ho fc Ohashi1ll998l : lTang et al.ll2009h . 



For the W51e2 region, we can identify two outflows: a high- velocity outflow and a 
low-velocity outflow that arise separately from the e2-E and e2-NW cores, respectively. 
We find that the northwest to southeast direction of the bipolar H 2 maser outflow from 
W51e2- E is in good agreement with the molecular outflow direction s een in the CO (2- 
1) line fceto fc Klaassen I l2008h and the CO(3-2) and HCN(4-3) lines Jshi. Zhao fc Han 



2010al lbh. but with higher expansion velocities for H2O masers. We classified maser features 
near W51e2 position into e2-E and e2-NW groups by their radial velocities: features with 
radial velocities within the rage of < Vlsr < 120 km s _1 were included in e2-NW and 
high-velocity features outside of the range in e 2-E. For comparison, t he m olecular radial 
velocity of W51e2 was found to be 55 km s _1 by IZhang. Ho fc Ohashi I (119981 ) from CH3CN 
(8-7) line observations. 



Discussion 



4.1. Distance to the Galactic Center 



Since the location of W51 is very close to the tangent point of the Sagittarius spiral arm 
of the Galaxy, the distance to W51 directly yields a good estimate of the distance R to the 
Galactic center by simple geometry. 

Figure |5] gives a schematic depiction of the source locations. At the galactic coordinates 
of (/, b) = (4949, —039), W51 lies very close to the tangent point (i.e., the angle 7 between 
the Sun and the Galactic center as viewed from W51 is very close to 90°; see Figure [5]), since 
the radial velocity of W51 Main/South, Vlsr ~ 58 km s -1 (Section l4.2p . is very close to the 
maximum radial velocity (which occurs at the tangent point) given by t> max = Oo(l — sinZ), 
assuming a flat rotati on curve of the Galaxy (i.e., d9/d-R = 0). For example, if we adopt 



R = 8.4 ± 0. 4 kpc flGhez et all 120081: iGillessen et al.l 120091) and Q /Ro = 29.45 ± 0.15 
km s^ 1 kpc -1 ( IReid fc Brunthaler 1120041 ) to obtain O , then the maximum velocity for I = 
4949 is f max = 59 ±3 km s , which is in good agreement with the observed radial velocity of 



W51 Main/South ((Vlsr) = 58 ±4 km s -1 ; see Section I4~2l) . Note that correction for a likely 
co unter-rotating pecu liar (non-circular) motion of a massive star forming region, as reported 
by IReid et al.l (j2009d ). would only make the observed radial velocity due to the Galactic 
rotation larger than the expected maximum velocity, which would then lend further support 
to the assumption that W51 lies at the tangent point. 
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From simple geometry, Rq = dcosb/ cos/, assuming only that W51 lies at the tangent 
point (i.e., 7 = 90°), our distance estimate of d = 5.41+°'^ kpc to W51 directly gives a 
trigonometric distance of Rq = 8.32^.0 43 ~ 8-32 ± 0.46 kpc to the Galactic center. 



We now relax the assumption that W51 is exactly at the tangent point and estimate 
the systematic effect this has on the uncertainty of Rq. For example, if we adopt Galactic 
and source parameters as discussed in Section 14.31 (i.e., 6 = 240 ± 20 km s _1 ; U@ = 
10.00 ± 0.36 km s _1 , V@ = 8.0 ± 3.0 km s" 1 , W & = 7.17 ± 0.38 km s _1 ; the source proper 
motion fi x = —2.64 ± 0.16 mas yr -1 , fi y = —5.11 ± 0.16 mas yr -1 ; the source radial velocity 
Vlsr = 58 ± 4 km s _1 and the r adial peculiar motion V src = — 4 ± 5 km s" 1 ) and 

solve for Rq using the equation by lReid. Moran fc Gwinn I ( 119881 ) with no assumption on the 
tangent point, we obtain R = 8.1 ± 1.1 kpc, which includes both statistical and systematic 
uncertainties. Therefore our best estimate of Rq from the parallax of W51 is Rq = 8.3 ±0.46 
(statistical) ±1.0 (systematic) kpc. Combining the statistical and systematic uncertainties in 
quadrature, we find Rq = 8.3 ±1.1 kpc. This estimate is consistent with other measurements 



of Rq: e.g., Rq = 8.0 ± 0.5 kpc by iReid I (jl993l ) from an en semble of cl a ssical techniques, 
and recent direct measurements of R n = 8.4 ± 0.4 kpc by iGhez et al.l (120081 ) and R - 
8.3 3 ±0.35 kpc by"" 

by 
by 



Reid et al. 



Reid et al. 



Gi 



lessen et al.l (120091 ) from stellar orbits around Sgr A*, R = 7.9 _ 7 kpc 



(2009b) from the trigonometric parallax of Sgr B2, and Rq = 8.4 ± 0.6 kpc 
(l2009ch from trigonometric parallaxes and proper motions of 16 massive star- 
forming regions across the Galaxy. We thus determined the distance to the Galactic center 
with good accuracy from a precise parallax measurement of W51 Main/South. 



4.2. Origins of Internal Maser Motions 

The internal proper motions of the H2O masers in W51 Main/South (Figure [2]) suggest 
four bipolar outflows, presumably from four young, massive stars: e2-NW, e2-E, e8, and u. 
We modeled the first three as uniformly expanding flows and fitted the data to obtain seven 
global parameters: the expansion speed and the three spatial and velocity components of the 
center of expansion. (There were not enough measured motions for the "u" group to reliably 
model.) For each maser feature, we fitted the measured three components of velocity using 
their positions on the sky as independent variables. Since the spatial offset along the line of 
sight from the center of expansion for each maser feature is unknown, we needed to solve for 
an additional offset parameter for each feature. 



We adopted a Bayesian fitting procedure using a Markov chain Monte Carlo method 
to explore parameter space. H2O m asers are known to have an internal dispersion of ~ 15 
km s _1 for each velocity component (IReid et al.lll9 88) and we added this in quadrature with 
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the (usually much smaller) measurement errors. Anticipating that some maser features may 
be "outliers" (i.e., not well modeled by simple uniform expansion), we adopted an "error- 
tolerant" probability density function for the data errors with the form probfoj crn) = Vo/v 2 , 
where a is the adopted (minimum) data uncertainty (see lSivia &: Skilling 1120071 ). For priors, 
we set the expansion speed (V exp ) at 75 ± 200 km s -1 ; the x and y proper motions of the 
center (Av x and Av y ) at ± 10 km s -1 relative to the reference feature #32; the LSR 
velocity of the center (Vlsr) was taken as 57 ± 10 km s _1 , a value between the average of 
th e thermal molecular lines in e 2 and e8 (Vlsr — 55 and 59 km s _1 for e2 and e8, respectively, 
by lZhang. Ho fc Ohashi Hl998l from the CH3CN (8-7) line); the coordinates of the center of 
expansion (x and y) were estimated visually from Figure |2] and assigned an uncertainty of 
±1 arcsec (i.e., (2, y) — (0 ± 1, 6.5 ± 1) arcsec for e2-NW, (0.3 ± 1, 6.2 ± 1) arcsec for e2-E, 
and (—0.5 ± 1, ± 1) arcsec for e8); and we assumed a source distance given by the parallax 
(7r) of 0.185 ± 0.010 mas determined in Section I3TT1 

The results of the fitting for each center of expansion are given in Table |5j We adopted 
the peak value of the histogram of each parameter as the best estimate of the parameter, 
and estimated its uncertainty from the halfwidth of the histogram where the number of 
occurrences becomes exp(— 1/2) of the peak, i.e., lcr for a Gaussian distribution. The error- 
weighted means of the fitted velocities are (Av x ) = —3 ±4 km s -1 and (Av y ) = 1 ±4 km s _1 
relative to feature 32, and (Vlsr) = 58 ± 4 km s -1 . These values are in good agreement 
with the unweighted mean of all 37 maser features in Table [3j (Av x ) = km s" 1 and 
(Av y ) = 2 km s^ 1 (at a distance of 5.41 kpc) and (Vlsr) = 62 km s _1 . We adopt the result 
from the Bayesian method as the best estimate for the systemic velocity of W51 Main/South. 
By adding (Av x ) = —3 ± 4 km s^ 1 and (Av y ) = 1 ± 4 km s" 1 (at a distance of 5.41 kpc) 
to the absolute proper motion of the feature 32, (n x , n y ) = (—2.53 ± 0.02, — 5.15 ± 0.04) 
mas yr _1 (see Section [331 and Table HJ), we obtain the systemic absolute proper motion of 
W51 Main/South to be (fi x ) = —2.64 ± 0.16 mas yr _1 and (/j, y ) = —5.11 ± 0.16 mas yr _1 , 
which will be used in Section 14.31 for fitting Galactic parameters. Figure [2] shows the maser 
motions relative to this systemic velocity. 



For comparison, Ku et al.l ( 120091 ) report proper motions of (fi x ,fi y ) = (—2.49 ± 0.07, 
-5.45 ± 0.14) mas yr" 1 and (-2.48 ± 0.08, -5.56 ± 0.08) mas yr" 1 for two CH 3 OH maser 
spot s in W51 IRS 2, whos e sky position is l'O west and 0'6 north of W51 Main/South (see, 
e.g., iFigueredo et al.ll2008l for a map of the giant H II region W51 A). At the source distance 
of 5.4 kpc, the two regions are separated by a sky-projected distance of ~ 1.8 pc. The 
eastward proper motions fi x of the two regions agree well, and the northward proper motions 
fi y differ by ~ 0.4 mas yr -1 ?^ 10 km s -1 . The difference in v y of 10 km s _1 is a reasonable 
value, considering the virial motion of ~ 7 km s -1 in each coord inate expected for each 
massive star that is exciting the maser emission ( iReid et al.l l2009d ) and some contribution 
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from the internal motion of the CH3OH maser spots « 3 km s" 1 (IMoscadelli et al.l 12002( 1 . 



Each H 2 maser outflow in W51 Main/South shows a relatively small range of inter- 
nal 3D speeds. This i s contrary to previous observations of other Galactic sta r-forming re- 



gions , e.g., Orion-KL ( iGenzel fc Downes Hl977tlGenzel et al.lll981al ) and W49N (IGwinn et al. 



19921 ). which suggested that high- and low- velocity outflows might originate from the same 



young stellar object. In W49, a strong acceleration region at a radius of 0.1 pc is required if 
all masers are associated with a single young stellar object. However, these conclusions were 
based on observations with less precise maser motions and, most importantly, much lower 
resolution dust emission maps. Our results for W51 suggest that multiple speed outflows 
may come from separate young stellar objects that can be very closely spaced on the sky. 



4.3. Galactic 3D Motion of W51 



Using our parallax distance and the fitted full-space systemic velocity of W51 Main/South 
obtained in Section |4T2| we now estimate the 3D motion of W51 Main/South in the Galaxy. 
As in Section I4.2[ we adopted a Bayesian fitting procedure using a Markov chain Monte 
Carlo method. We set priors for nine parameters: the source distance of W51 Main/South 
given by the parallax (ir) determined in Section 13. 1( distance to the Galactic center (i? ) 
and the angular rotation speed (Qq/Rq) of the Galaxy at the LSR; the Solar motion (U Q , 
Vq, W & ) relative to the LSR; and source peculiar motion (relative to a circular orbit) of 
W51 Main/South (U SIC , V 8TC , W SIC ), where U, V, and W denote velocity components in the 
directions (at the source location) toward the Galactic center, toward Galactic rotation and 
toward the north Galactic pole (NGP), respectively. 



The prior R = 8.25 ±0.4 kpc was chosen a s an average between 8.0 kpc bvlReid I (119931 ) 



and the recent direct measurements of Rq (e.g.. 
prior value of 0q/-Ro = 29.45±0.15 k m s" 1 kpc" 1 (for Vjp 



Ghez et al.ll2008l : iGillessen et al.ll2009h . The 



5.25 km s 1 and corrected for the 



varied V® values) was adopted from iReid fc Brunthaler I (120041 ) from the proper motion of 
Sgr A*. The priors of U Q = 10.00 ±0.36 km s " 1 ari d W & = 7.17 ± 0.38 km s" 1 were adopted 
from Hipparcos values by iDehnen fc Binnev I (119981) . while V@ = 8.0 ± 3.0 km s _1 was chosen 
to allow for both 5. 25 km s~ l bvlDehnen fc Binnev I (119981 ) and a larger value of ~ 11 km s _1 
recently revised by lBinney I ( 120101 ). We set the priors of U STC , V STC and W src to be rather open, 



± 20 km s \ in particu l ar to allow for a wide possible range of V^ c due to the uncertaint y 
of V Q jReid et al.l l2Q09cl : iBinnev I bold : iMcMillan fc Binnev I boiol ; ISchonrich etHlboioh . 



We assumed a flat rotation curve of the Galaxy (i.e., dQ/dR = 0). The nine parameters 
were then adjusted to maximize the probability of the model given the above priors and 
the data of the source proper motion and LSR velocity, (fi x ) = —2.64 ± 0.16 mas yr _1 , 
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(fj, y ) = —5.11 ± 0.16 mas yr _1 and (Vlsr) = 58 ± 4 km s _1 . Note that the proper motions 
(fi x , fly) are in the heliocentric frame, while the given radial velocity Vlsr is the LSR velocity, 
converted from the heliocentric radial velocity fheiio (obtained fr om the Doppler sh ift) using 
the standard solar motion by definition relative to the LSR (see lReid et al.ll2009q ). 



Since the obtained probability distribution functions for some parameters did not follow 
a Gaussian function, we adopted, as a conservative approach, the mean value of each pa- 
rameter as the best estimate for the parameter and the standard deviations from the mean 
as the uncertainty, based on the last 90% of 500,000 trials. Most of the Galactic parame- 
ters returned the uncertainties only slightly smaller than the priors, indicating that these 
parameters were mainly constrained by the fairly strong priors. 

The fitted result yields (U SIC , V SIC , W SIC )= (-7 ±5, -4 ± 5, 4 ± 4) km s" 1 for the 
peculiar (non-circular) mot ion of W5 1 Mai n/South. For comparison, the parallax and full- 
space motion reported by Ku et al.l ( 120061 ) for W3(OH) in the Perseus spiral arm would 
yield (U STC , V STC , W STC )= (19 ± 3, —15 ±4, 1 ± 2) km s _1 under the same priors for Galactic 
parameters. Therefore, the full-space velocity of W51 Main/South suggests that W51 is in 
a nearly circular orbit about the Galactic center, with no large peculiar motion. 
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Fig. 1. — 22 GHz H2O maser spectrum of W51 Main/South over all four 8 MHz frequency 
bands from the first-epoch data. Upper panel: on full scale. Lower panel: blowup for 
weaker emission. Dotted vertical lines indicate the boundaries of four 8 MHz frequency 
bands observed. 
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Fig. 2.— Left panel: map of the positions and internal motions of H2O maser fea- 
tures in W51 Main/South, as seen from the mean motion of all regions (see Section 
142]) . Absolute coordinates of the map origin are R.A. (J2000)= 19 h 23 m 43?93427 and decl. 
(J2000)= 14°30 / 28('3498 (i.e., position of reference feature 32 at the mid-observation year of 
2009.318). The HC and UC H II regions are marked by crosses with positions and beam sizes 
from previous observations (see the text for references). A sample motion of 2.0 mas yr -1 
(«51 km s -1 at a distance of 5.4 kpc) is shown at the bottom left corner of the map. Right 
panel (top): blowup for W51e2-NW region with lower- velocity maser spots. Right panel 
(bottom): H 2 maser spectra of the four regions with intensities summed over all spots 
detected at the first epoch. 
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Fig. 3. — Images of background continuum sources. Source names are shown in the upper left 
corner and restoring beams are indicated by gray ellipse in the lower left corner of each panel. 
The more heavily resolved source, J1922+1504, was not used for parallax measurement of 
W51 Main/South. The contours are at 10%, 20%, 30%, 90% of the peak intensities, 
which are 0.108 Jy beam" 1 for J1922+1530, 0.633 Jy beam" 1 for J1924+1540, and 0.025 
Jy beam" 1 for J1922+1504. 



Table 1. VLB A Observations of W51 Main/South 



Epoch 


Date 


Day of Year 


Time Range 


Restoring Beam 


Antennas Unavailable 








(UT) 


(mas, mas, deg) 


No Data 


Flagged 


1 


2008 Oct 22 


2008/296 


19:36-04:31 


1.2x0.4 at -20 


SC 


HN 


2 


2009 Apr 27 


2009/117 


07:20-16:16 


0.8x0.3 at -7 




NL 


3 


2009 Apr 30 


2009/120 


07:09-16:04 


1.1x0.3 at -18 


KP 


NL, SC 


4 


2009 Oct 20 


2009/293 


19:44-04:40 


1.1x0.3 at -19 


SC 





Note. — Restoring beam sizes (FWHM) are indicated in Column 5 by the major and minor 
axes in mas and position angles (P.A.) east of north in degrees. Antenna codes are HN: Hancock, 
KP: Kitt Peak, NL: North Liberty, and SC: Saint Croix. 
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Fig. 4. — Parallax and proper motion data (points with error bars) and fits (solid or dashed 
curves) for W51 Main/South. Plotted are position measurements of an H 2 maser spot at 
Vlsr = 62.9 km s -1 relative to background quasars J1922+1530 and J1924+1540 with first 
and fourth epochs labeled. Left panel: sky-projected motion of the maser spot relative to 
the two quasars. Note that the second and third epochs were spaced only by 3 days and 
thus indistinguishable to eye in the plot. Middle panel: east and north position offsets and 
the best-fit parallax and proper motion of the maser spot vs. time. Solid and dashed curves 
indicate fits for eastward and northward positions, respectively. Constant positional offsets 
from data are added in each coordinate for clarity. Right panel: same as the middle panel, 
except the best-fit proper motion has been removed to show the effects of the parallax only. 
Filled triangles show maser positions measured relative to J1922+1530 and open circles to 
J1924+1540. The northward data have been offset from the eastward data for clarity. 
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Fig. 5. — Schematic depiction of source locations and Galactic parameters when the source 
is at the "tangent point" (7 = 90°) in the Galaxy. 
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Table 2. Source Positions and Brightnesses 



Source 


R.A. (J2000) 


Dec. (J2000) 


a 

"sep 


P.A. 


Tb Vlsr 




(h m s) 


(d ' ") 


(deg) 


(deg) 


(Jy/b) (km s" 1 ) 


W51 Main/South 


19 23 43.87363 


+14 30 29.4536 






746 58 


J1922+1530 


19 22 34.699314 


+ 15 30 10.03262 


1.0 


-16 


0.108 


J1924+1540 


19 24 39.455870 


+ 15 40 43.94172 


1.2 


11 


0.633 


J1922+1504 


19 22 33.27291 


+15 04 47.5382 


0.6 


-26 


0.025 



Note. - Coordinates are those used in the VLBA correlator. Absolute positions 
of J1922+1530 an d J1924+1540 are based on the VLBA Calibrator Survey 1 (VCS1; 



Beasley et al.l 120021 ) and known in each coordinate to about ±0.3 mas (J1922+1530) and 



±0.8 mas (J1924+1540). The position of J1922+1504 is accurate to about ±20 mas. Angu- 
and position angles (P.A.) east of north relative to the H 2 maser source in 



'sep, 



lar offsets 

W51 Main/South are indicated in Columns 4 and 5. Peak brightnesses (T& 
and background continuum sources are from the first epoch. 



for the maser 
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Table 3. Maser Feature Motions in W51 Main/ South 



ID 






X 


yo 


H 


X 






(km s J 


(Jy/b) 


(' 

\ 


') 


(' 


) 


(mas yr x ) 


(mas yr -1 ) 


^ez-iN vv j 






















i 
i 


oy .o 


1 .0 


-0.879 


(0.008) 


7.786 


(0.025) 


-0.43 


(0.03) 


1.31 


(0.11) 


o 
L 


4^ n 


ZL 9 


0.721 


(0.008) 


6.177 


(0.008) 


0.14 


0.02 

V / 


-0.23 


0.02 

V / 


Q 
O 


4fi fi 
^u.u 


Q 3 
y .o 


0.700 


(0.005) 


6.196 


(0.003) 


0.75 


(0.01) 


-0.49 


(0.01) 


4 


^1 4 


fi n 


0.930 


(0.023) 


7.216 


(0.001) 


-0.32 


(0.06) 


-0.14 


(0.01) 


c; 

(J 


fin 7 


fi44 ^ 


-0.732 


(0.010) 


7.238 


(0.005) 


-0.60 


(0.02) 


-1.15 


(0.01) 


u 




Q 

oou.y 


-0.732 


(0.006) 


7.238 


(0.024) 


0.06 


(0.02) 


-0.87 


(0.06) 


7 




1 14 Q 


-0.674 


(0.003) 


7.518 


(0.003) 


-0.53 


(0.01) 


-0.04 


(0.01) 


o 


77 ^ 


9 ^ 


-0.674 


(0.044) 


7.315 


(0.023) 


-0.85 


(0.11) 


-0.14 


(0.05) 


q 


oo.u 




-0.662 


(0.020) 


7.416 


(0.021) 


-0.41 


(0.09) 


0.14 


(0.10) 


1 D 




^7 ^ 


-0.655 


(0.003) 


7.367 


(0.002) 


-0.89 


(0.01) 


1.37 


(0.01) 


1 1 


i nn 1 


HOI 3 


-0.656 


(0.006) 


7.361 


(0.003) 


-0.21 


(0.01) 


1.30 


(0.01) 


1 9 


1UU.O 


7^ zL 


-0.702 


(0.007) 


7.306 


(0.012) 


-0.51 


(0.04) 


1.11 


(0.12) 


1 3 


1 HQ ^ 


1 9 


-0.701 


(0.061) 


7.307 


(0.087) 


1.00 


(0.13) 


3.11 


(0.18) 
























14 


-51.7 


3.6 


0.800 


(0.011) 


6.033 


(0.030) 


5.64 


(0.02) 


-3.74 


(0.07) 


1 u 


— ^ Q 


18 1 

lO. 1 


0.706 


(0.013) 


6.059 


(0.006) 


5.01 


(0.04) 


-2.90 


(0.02) 


1 fi 


—93 fi 


34 9 


0.687 


(0.003) 


6.095 


(0.005) 


5.20 


(0.02) 


-2.98 


(0.01) 


1 7 


1Q7 n 

lO ( .y 


4 ^ 

^.O 


-0.210 


(0.015) 


6.683 


(0.003) 


-3.22 


(0.04) 


2.26 


(0.01) 


lo 


100.0 


9 7 


0.269 


(0.008) 


6.323 


(0.011) 


-5.24 


(0.02) 


1.42 


(0.02) 


1 Q 
i y 


1 3Q fi 


1 n 

1 .u 


-0.229 


(0.008) 


6.670 


(0.004) 


-3.53 


(0.03) 


2.18 


(0.01) 


20 


139.8 


12.0 


0.269 


(0.016) 


6.323 


(0.013) 


-4.65 


(0.04) 


2.01 


(0.03) 


21 


149.8 


6.2 


-0.646 


(0.007) 


7.125 


(0.010) 


-2.83 


(0.03) 


2.98 


(0.03) 


22 


150.9 


3.5 


-0.649 


(0.017) 


7.083 


(0.017) 


-2.64 


(0.05) 


2.61 


(0.05) 


(e8) 

23 


-28.7 


10.0 


0.591 


(0.001) 


-0.311 


(0.005) 


6.78 


(0.01) 


-3.49 


(0.02) 


24 


-7.4 


5.5 


-0.938 


(0.002) 


0.034 


(0.010) 


-3.90 


(0.01) 


2.27 


(0.02) 


25 


26.6 


6.7 


-0.198 


(0.005) 


0.049 


(0.007) 


1.84 


(0.02) 


0.25 


(0.02) 


26 


41.3 


81.7 


0.365 


(0.073) 


-0.074 


(0.075) 


3.01 


(0.15) 


-0.80 


(0.16) 


27 


43.0 


37.3 


-1.404 


(0.003) 


0.039 


(0.005) 


-2.81 


(0.01) 


0.66 


(0.02) 
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Table 3 — Continued 



ID 


Vlsr 


T 


Xq 






Ms/ 




(km s -1 ) 


(Jy/b) 


(") 


(") 


(mas yr -1 ) 


(mas yr -1 ) 


28 


42.3 


11.5 


0.365 (0.027) 


-0.074 (0.020) 


1.93 (0.07) 


0.13 (0.06) 


29 


52.1 


27.8 


0.172 (0.012) 


-0.040 (0.023) 


1.52 (0.03) 


0.31 (0.06) 


30 


60.0 


63.6 


0.277 (0.003) 


-0.018 (0.001) 


0.40 (0.01) 


0.07 (0.01) 


31 


61.3 


128.0 


-0.534 (0.003) 


0.534 (0.006) 


0.50 (0.02) 


-0.31 (0.02) 


32 


62.9 


918.1 


0.000 (0.001) 


0.000 (0.001) 


0.00 (0.01) 


0.00 (0.01) 


33 


63.2 


71.9 


-0.905 (0.002) 


-1.579 (0.004) 


0.35 (0.01) 


-0.19 (0.01) 


(u) 
34 


-25.7 


78.5 


-2.054 (0.105) 


-8.488 (0.119) 


-0.14 (0.22) 


-0.27 (0.25) 


35 


40.3 


15.3 


-2.767 (0.011) 


-7.023 (0.010) 


-0.16 (0.03) 


-3.02 (0.03) 


36 


42.7 


1.7 


-2.765 (0.015) 


-7.022 (0.010) 


-0.11 (0.04) 


-2.88 (0.03) 


37 


111.9 


107.2 


-2.787 (0.019) 


-6.876 (0.023) 


-0.78 (0.05) 


1.62 (0.05) 



Note. - - H 2 maser features in W51 Main/South listed with relative positions x , 
yo and proper motions /i x , \i y with respect to feature 32. x ad y denote eastward and 
northward directions, respectively. A proper motion of 1 mas yr _1 corresponds to 25.6 
km s _1 at a distance of 5.41 kpc. The uncertainties were obtained from the formal fitting 
uncertainties and are shown in parentheses. The positional offsets of x , y relative to 
feature 32 are those fitted for the mid-observation year of 2009.318 (i.e., average year 
of all four epochs). The brightness (TJ,) of each maser feature gives the maximum value 
of all spots of the feature from all four epochs. 
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Table 4. Parallax and Proper Motion Fits for W51 Main/South 



Maser Vlsr 


Background 


Parallax 


fix 




(km s -1 ) 


Source 


(mas) 


(mas yr -1 ) 


(mas yr -1 ) 


62.9 


J1922+1530 


0.183 ±0.006 


-2.56 ±0.02 


-5.09 ±0.03 


62.9 


J1924+1540 


0.187 ±0.009 


-2.49 ±0.02 


-5.21 ±0.04 


62.9 


Combined 


0.185 ±0.007 


-2.53 ±0.02 


-5.15 ±0.04 



Note. - Combined fit used a single parallax parameter and a single 
proper motion in each coordinate for the maser spot position relative to 
the two background sources. The uncertainty for each parallax fit was 
obtained from the formal fitting uncertainty. 



Table 5. Fitted Motions of Proto-stars and Expansion Speed 



Source 


X 


y 


Av x 


AVy 




^exp 




(") 


(") 


(km s _1 ) 


(km s _1 ) 


(km s _1 ) 


(km s _1 ) 


e2-NW 


-0.7 ±0.6 


6.3 ±0.8 


-8 ±6 


2±5 


66 ±7 


20 ± 9 


e2-E 


2.0 ±0.4 


5.8 ±0.5 


-2 ±6 


1 ±9 


60 ±8 


120 ± 12 


e8 


-1.1 ±0.6 


-0.1 ±0.7 


6±8 


1 ±6 


49 ±7 


21 ± 11 



Note. - 



(x, y) and (Av x , Av y ) are relative to the reference feature 32. 
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Table 6. Maser Spot Motions in W51 Main/South 



ID 


# 


^LSR 

(km s _1 ) 


(Jy/b) 


x 
(") 


(") 


(mas yr _1 ) 


(mas yr _1 ) 


1 


1 


on n 

oy.U 


i.i 


-0.879 


(0.015) 


7.786 


(0.046) 


-0.41 


(0.05) 


1.35 


(0.14) 


1 


o 
Z 


on a 

oy.4 


l.o 


-0.879 


(0.010) 


7.786 


(0.030) 


-0.45 


(0.05) 


1.23 


(0.17) 


Z 


1 
1 


A A Q 
44. 


A 1 
4.Z 


0.721 


(0.009) 


6.177 


(0.008) 


0.15 


(0.02) 


-0.25 


(0.02) 


Z 


Z 


40. o 


Z. ( 


0.721 


(0.026) 


6.177 


(0.025) 


-0.02 


(0.06) 


-0.10 


(0.05) 


9 

o 


1 

1 


4o.o 


O Q 

y.o 


0.700 


(0.015) 


6.196 


(0.003) 


0.59 


(0.03) 


-0.49 


(0.01) 


Q 

o 


o 
z 


/l P Cl 

4o.y 


K 7 
0. ( 


700 

\J • I \J\J 


(0 0061 


6 196 


(0 0071 


78 


fO 01 1 


-0 49 


(0 Oil 


4 


-1 

1 


01. z 


k p 
0.0 


0.930 


(0.031) 


7.216 


(0.012) 


-0.34 


(0.08) 


-0.17 


(0.03) 


4 


o 
z 


£1 p 
01.0 


O.U 


0.930 


(0.036) 


7.216 


(0.000) 


-0.29 


(0.10) 


-0.14 


(0.01) 





-1 

1 


7 

DO. 1 


1U0.O 


-0.732 


(0.040) 


7.238 


(0.031) 


-0.67 


(0.11) 


-1.25 


(0.08) 





z 


oy.z 


1 07 7 
Iz 1 . 1 


-0 732 


(0 0381 


7 238 


(0 0271 


-0 71 


fO 091 


-1.25 


(0 071 


r 




Q 

o 


Oy.O 


1 t^P Q 
100. 


-0 732 


(0 0251 


7 238 


(0.022) 


-0 59 


fO 051 


-1 16 


(0 051 





4 


OU.U 


olO.O 


-0.732 


(0.016) 


7.238 


(0.007) 


-0.53 


(0.04) 


-1.11 


(0.02) 


r 



r 



oU.4 


ro7 £ 
Oo 1 .0 


-0 732 


(0 0201 


7 238 


(0 0071 


-0 67 


fO 051 


-1 18 


(0 021 








PC\ Q 

OU.o 


f; /i /i k 
044.0 


-0.732 


(0.052) 


7.238 


(0.031) 


-0.68 


(0.12) 


-1.16 


(0.07) 







P 1 9 
Ol.O 


/I 7/1 7 

4(4. ( 


-0.732 


(0.048) 


7.238 


(0.040) 


-0.44 


(0.21) 


-1.11 


(0.13) 


c 



1 


£1 7 
01. ( 


ooO.y 


-0.732 


(0.006) 


7.238 


(0.039) 


0.02 


(0.03) 


-0.94 


(0.18) 





z 


PO 1 

0Z.1 


zzi.y 


-0.732 


(0.016) 


7.238 


(0.040) 


0.12 


(0.04) 


-0.87 


(0.09) 





Q 

o 


Oz.O 


loo.O 


-0.732 


(0.027) 


7.238 


(0.048) 


0.06 


(0.06) 


-0.83 


(0.10) 


17 
I 


1 

1 


(l.o 


70 Q 

1 y.o 


-0.674 


(0.006) 


7.518 


(0.013) 


-0.58 


(0.02) 


0.04 


(0.08) 


iy 
t 


o 

z 


71 1 
( Z.Z 


11A 

11U.0 


-0.674 


(0.006) 


7.518 


(0.011) 


-0.54 


(0.01) 


-0.01 


(0.04) 


n 
t 


Q 

o 


71 P 

1 Z.O 


1 1 /i n 
114. y 


-0.674 


(0.006) 


7.518 


(0.007) 


-0.53 


(0.03) 


-0.02 


(0.02) 


7 


4 


73.1 


91.1 


-0.674 


(0.007) 


7.518 


(0.008) 


-0.51 


(0.04) 


-0.03 


(0.03) 


7 


5 


73.5 


48.5 


-0.674 


(0.006) 


7.518 


(0.008) 


-0.49 


(0.02) 


-0.07 


(0.03) 


7 


6 


73.9 


17.8 


-0.674 


(0.013) 


7.518 


(0.007) 


-0.52 


(0.03) 


-0.06 


(0.02) 


8 


1 


77.3 


2.1 


-0.674 


(0.067) 


7.315 


(0.036) 


-0.66 


(0.18) 


-0.07 


(0.09) 


8 


2 


77.7 


2.5 


-0.674 


(0.059) 


7.315 


(0.030) 


-0.94 


(0.13) 


-0.17 


(0.07) 


9 


1 


82.7 


6.8 


-0.662 


(0.035) 


7.416 


(0.054) 


-0.39 


(0.22) 


0.16 


(0.31) 


9 


2 


83.2 


8.8 


-0.662 


(0.043) 


7.416 


(0.041) 


-0.31 


(0.23) 


0.16 


(0.21) 


9 


3 


83.6 


13.5 


-0.662 


(0.048) 


7.416 


(0.040) 


-0.44 


(0.16) 


0.14 


(0.15) 


9 


4 


84.0 


13.6 


-0.662 


(0.034) 


7.416 


(0.035) 


-0.45 


(0.14) 


0.11 


(0.17) 
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Table 6 — Continued 



ID 


# 


Vlsr 

n — 1 \ 
(km s J 


T b 

/ T /I \ 

(Jy/b) 


x 
(") 


(") 


(J>x 

(mas yr _1 ) 


(mas yr _1 ) 


ID 


i 

X 


Q4 Fi 


93 8 

<_/W . (J 


-n.655 


(0.010) 


7.367 


(0.018) 


-0.90 


(0.03) 


1.43 


(0.05) 


1 n 


9 

Zj 


Q4 Q 


2n^ i 

ZjWW. X 


-n.655 


(0.010) 


7.367 


(0.012) 


-0.88 


(0.04) 


1.40 


(0.05) 


1 n 


(J 


4 

C/W .t: 


^7 1 


-n.655 


(0.012) 


7.367 


(0.008) 


-0.86 


(0.03) 


1.41 


(0.03) 


1 n 


4 


<vW . O 


3^)7 Fi 


-n.655 


(0.005) 


7.367 


(0.005) 


-0.85 


(0.01) 


1.41 


(0.01) 


1 n 


W 


Q6 2 

<-/W . Zj 


1 82 4 

X O Zj . t: 


-n.656 


(0.006) 


7.367 


(0.003) 


-0.98 


(0.02) 


1.36 


(0.01) 


1 n 


6 

W 


Q6 6 

<-/W . W 


51 7 

W X . 1 


-n.656 


(0.010) 


7.367 


(0.010) 


-1.06 


(0.03) 


1.26 


(0.08) 


1 1 


1 

X 


i nn n 


UU X . o 


-n.656 


(0.024) 


7.361 


(0.005) 


-0.18 


(0.05) 


1.28 


(0.01) 


X X 


9 

Zj 


i nn 4 

X WW . t: 


28n 6 

Zj(J w . w 


-n.656 


(0.021) 


7.361 


(0.005) 


-0.17 


(0.05) 


1.30 


(0.01) 


1 1 


W 


i nn 8 

X WW . o 


76 Q 


-n.656 


(0.010) 


7.361 


(0.010) 


-0.21 


(0.02) 


1.33 


(0.02) 


1 1 


4 

t: 


im 3 


2n 1 

Zj W . X 


-n.656 


(0.008) 


7.361 


(0.014) 


-0.22 


(0.02) 


1.36 


(0.03) 


1 2 


1 

X 


i n^ q 

X WW . u 


2Q 7 


-n.7n2 


(0.024) 


7.306 


(0.044) 


-0.56 


(0.11) 


1.22 


( n -21) 


12 


2 


m6 3 

X WW . KJ 


48 4 

t:0 . t: 


-n.7n2 


(0.036) 


7.306 


(0.034) 


-0.52 


(0.12) 


1.12 


(0.26) 


12 

XZj 




1 n6 7 

X WW . 1 


63 n 

WO . w 


-n.7n2 


(0.017) 


7.306 


(0.031) 


-0.43 


(0.13) 


1.00 


(0.28) 


1 2 


4 


i n7 2 

XW I .Zj 


73 4 

I KJ .t: 


-n.7n2 


(0.009) 


7.306 


(0.016) 


-0.50 


(0.05) 


1.07 


(0.20) 


1 3 


1 

X 


1 ns 8 


^2 1 

WZj . X 


-0.701 


(0.082) 


7.307 


(0.134) 


1.03 


(0.17) 


3.17 


(0.27) 


1 3 

X w 


2 

Zj 


1 nn 7 


in3 2 

X W W . Zj 


-0.701 


(0.092) 


7.307 


(0.114) 


0.97 


(0.19) 


3.07 


(0.24) 


14 

It: 


1 

X 


-52 2 

WZj . Zj 


1 Q 

X . u 


o.snn 


(0.017) 


6.033 


(0.044) 


5.64 


(0.04) 


-3.74 


(0.09) 


14 

It: 


2 

Zj 


-51 8 

W X . (J 


3 6 

w . w 


n.snn 


(0.019) 


6.033 


(0.061) 


5.65 


(0.04) 


-3.74 


(0.14) 


14 

It: 


(J 


-51 4 

W X .t: 


3 3 

w . w 


n.snn 


(0.026) 


6.033 


(0.056) 


5.63 


(0.06) 


-3.75 


(0.12) 


1 Fi 

X W 


1 

X 


-35 4 

Ua . t: 


1 1 

X . X 


n.7n6 


(0.016) 


6.059 


(0.017) 


5.06 


(0.06) 


-2.83 


(0.06) 




2 

Zj 


-34 9 


2.1 


n.7n6 


(0.043) 


6.059 


(0.010) 


4.97 


(0.12) 


-2.87 


(0.03) 


15 


3 


-34.5 


6.1 


n.7n6 


(0.041) 


6.059 


(0.014) 


4.97 


(0.10) 


-2.95 


(0.04) 


15 


4 


-34.1 


16.8 


n.7n6 


(0.047) 


6.059 


(0.021) 


5.05 


(0.13) 


-2.95 


(0.05) 


15 


5 


-33.7 


18.1 


n.7n6 


(0.046) 


6.059 


(0.019) 


4.92 


(0.13) 


-2.90 


(0.05) 


16 


1 


-24.8 


in.i 


n.687 


(0.010) 


6.095 


(0.016) 


5.38 


(0.08) 


-3.40 


(0.10) 


16 


2 


-24.4 


22.3 


n.687 


(0.004) 


6.095 


(0.053) 


5.31 


(0.04) 


-3.26 


(0.35) 


16 


3 


-24.0 


23.2 


n.687 


(0.015) 


6.095 


(0.077) 


5.23 


(0.05) 


-2.98 


(0.33) 


16 


4 


-23.6 


3n.9 


n.687 


(0.010) 


6.095 


(0.005) 


5.14 


(0.03) 


-2.97 


(0.01) 


16 


5 


-23.2 


34.2 


n.687 


(0.018) 


6.095 


(0.070) 


5.07 


(0.05) 


-3.51 


(0.17) 


16 


6 


-22.7 


16.6 


n.687 


(0.024) 


6.095 


(0.071) 


5.06 


(0.13) 


-3.42 


(0.28) 
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Table 6 — Continued 



ID 


# 


(km s J 


T b 

( T /I \ 

(Jy/b) 


x 
(") 


(") 


fJ>x 

(mas yr _1 ) 


(mas yr 


) 


1 7 

X i 


1 

X 


1 37 7 


3 9 


-0.210 


(0.036) 


6.683 


(0.003) 


-3.22 


(0.09) 


2.26 


(o 


01 ) 

VI) 


1 7 

X i 


9 

Zj 


1 38 1 

X W(J . X 


4 ^ 

t: . W 


-0.210 


(0.016) 


6.684 


(0.015) 


-3.22 


(0.05) 


2.25 


(o 


94) 

V<i) 


1 8 


1 


1 38 1 


1 9 

X . Zj 


0.269 


(0.014) 


6.323 


(0.031) 


-5.23 


(0.05) 


1.45 


(o 


08) 

WO 1 


1 8 

X O 


2 

Zj 


1 38 5 

1U(J . W 


2 7 

Zj . 1 


0.269 


(0.010) 


6.323 


(0.011) 


-5.25 


(0.02) 


1.41 


(o 


02") 

WZj J 


1 9 


1 

X 


1 39 4 


1 9 

X . W 


-0.229 


(0.016) 


6.670 


(0.006) 


-3.50 


(0.05) 


2.17 


(o 


02) 

WZj J 


1 9 


2 

Zj' 


1 39 8 


9 9 


-0.229 


(0.009) 


6.670 


(0.007) 


-3.54 


(0.03) 


2.21 


(o 


02") 

WZj J 


90 

Zj W 


1 

J. 


1 3Q 4 


C/ . W 


0.269 


(0.048) 


6.323 


(0.035) 


-4.76 


(0.10) 


1.98 


(o 


07) 


20 

Zj w 


2 

Zj 


1 39 8 


1 2 9 

X Zj . W 


0.269 


(0.028) 


6.323 


(0.021) 


-4.67 


(0.06) 


2.00 


(o 


0^ 

WW 1 


29 

Zj w 


W 


149 2 

X t: W . Zj 


9 4 


0.269 


(0.031) 


6.323 


(0.022) 


-4.63 


(0.08) 


2.03 


(o 


WW 1 


29 

Zj w 


4 


149 6 

1t:U . W 


3 7 

W . 1 


0.269 


(0.028) 


6.323 


(0.049) 


-4.62 


(0.06) 


2.06 


(o 


11"! 


21 

Z. X 


1 


149 ^ 

1t:J . W 


W . W 


-0.646 


(0.014) 


7.125 


(0.011) 


-2.78 


(0.04) 


3.00 


(o 


04) 

V<i) 


21 


2 


149 9 


6 2 

V . Zj 


-0.646 


(0.010) 


7.125 


(0.026) 


-2.83 


(0.05) 


2.90 


(o 


15) 

lO) 


21 

Z X 




1 59 3 

X WW . w 


3 2 

W . Zj 


-0.646 


(0.017) 


7.125 


(0.048) 


-2.92 


(0.06) 


2.73 


(0 


1 6") 

IV) 


22 


i 

X 


1 59 7 

X WW . 1 


3 ^ 

W . W 


-0.649 


(0.036) 


7.083 


(0.037) 


-2.66 


(0.09) 


2.73 


(o 


10) 


22 

_ Z, 


2 

Zj 


1 ^1 2 

X W X . Zj 


2 3 

Zj . W 


-0.649 


(0.020) 


7.083 


(0.018) 


-2.62 


(0.06) 


2.57 


(o 


06) 

WW 1 


23 

Zj W 


1 

X 


-29 1 

Zj<_/ . X 


8 1 

(J . X 


0.591 


(0.003) 


-0.311 


(0.008) 


6.77 


(0.01) 


-3.50 


(o 


02) 

WZj J 


23 

Zj W 


2 

Z 


-28 6 

Zj(J . W 


19 9 


0.591 


(0.003) 


-0.311 


(0.009) 


6.78 


(0.01) 


-3.50 


(o 


03) 

WW J 


23 

Zj W 


w 


-28 2 

Zj(J . Zj 


^ 6 

W . W 


0.591 


(0.002) 


-0.311 


(0.013) 


6.78 


(0.01) 


-3.48 


(o 


04) 

V<i) 


23 

Zj W 


4 


-27 8 

Zj 1 .(J 


1 6 

X . W 


0.591 


(0.010) 


-0.311 


(0.019) 


6.79 


(0.04) 


-3.47 


(o 


07) 

VI ) 


24 

Zjt: 


1 

X 


-8 9 

O . W 


3 8 

W . O 


-0.938 


(0.004) 


0.034 


(0.025) 


-3.82 


(0.01) 


2.27 


(o 


07) 
u i ; 


24 


2 


-7 6 


5 5 

'-J . W 


-0.938 


(0.007) 


0.034 


(0.017) 


-3.85 


(0.02) 


2.19 


(o 


07) 

V 1 ) 


24 


3 


-7.2 


4.3 


-0.938 


(0.003) 


0.034 


(0.018) 


-3.93 


(0.01) 


2.20 


(o 


04) 


24 


4 


-6.3 


2.7 


-0.938 


(0.032) 


0.034 


(0.019) 


-3.85 


(0.07) 


2.34 


(o 


04) 


25 


1 


26.3 


4.5 


-0.198 


(0.008) 


0.049 


(0.009) 


1.84 


(0.02) 


0.25 


(o 


02) 


25 


2 


26.7 


6.7 


-0.198 


(0.006) 


0.049 


(0.012) 


1.82 


(0.05) 


0.23 


(o 


07) 


26 


1 


41.1 


54.6 


0.365 


(0.118) 


-0.074 


(0.099) 


3.00 


(0.26) 


-0.83 


(o 


21) 


26 


2 


41.5 


81.7 


0.365 


(0.094) 


-0.074 


(0.115) 


3.01 


(0.19) 


-0.76 


(o 


24) 


27 


1 


42.3 


16.3 


-1.404 


(0.003) 


0.039 


(0.017) 


-2.81 


(0.01) 


0.67 


(o 


04) 


27 


2 


42.7 


37.3 


-1.404 


(0.005) 


0.039 


(0.023) 


-2.83 


(0.01) 


0.64 


(o 


06) 


27 


3 


43.2 


36.9 


-1.404 


(0.018) 


0.039 


(0.022) 


-2.85 


(0.05) 


0.62 


(o 


06) 



-30 - 



Table 6 — Continued 



ID 


# 


(km s J 


T b 

( T /I \ 

(Jy/b) 


x 
(") 


(") 


fJ>x 

(mas yr 


) 


(mas yr -1 ) 


27 


4 


43 6 

TI-J . W 


16 8 

x w . o 


-1.404 


(0.044) 


0.039 


(0.014) 


-2.86 


(o 


1 2) 


0.64 


(0.07) 


27 


(J 


44 n 

11 . w 


15 5 

X W . W 


-1.404 


(0.029) 


0.039 


(0.009) 


-2.82 


(o 


19") 


0.66 


(0.04) 


97 


U 


44 4 

tit: .t: 


8 1 

O. J. 


-1.404 


(0.023) 


0.039 


(0.009) 


-2.90 


(o 


1 31 


0.68 


(0.05) 


28 

Zj O 


i 

X 


42 3 


1 1 ^ 

J. J. . w 


0.365 


(0.041) 


-0.074 


(0.031) 


1.98 


(o 


11) 


0.16 


(0.09) 


28 

Zj(J 


9 

Zj 


42 7 

t:Zj . 1 


2 4 

Zj . t: 


0.365 


(0.035) 


-0.074 


(0.027) 


1.90 


(o 


W<J 1 


0.11 


(0.07) 


2Q 

Zj U 


1 

X 


^2 n 

WZj . w 


27 8 

Zj 1 .(J 


0.172 


(0.018) 


-0.040 


(0.029) 


1.48 


(o 


041 

V^) 


0.34 


(0.07) 


2Q 

ZjC/ 


o 

Zj 


^2 4 

WZj . t: 


x w . w 


0.172 


(0.018) 


-0.040 


(0.039) 


1.57 


(o 


041 

U4: 1 


0.26 


(0.09) 


30 

o w 


1 

X 


5Q 6 

W<J . W 


37 7 

W I . I 


0.277 


(0.005) 


-0.018 


(0.004) 


0.31 


(o 


WZj J 


0.13 


(0.01) 


30 


9 

Zj 


6D n 

WW . W 


63 6 

WW . W 


0.277 


(0.007) 


-0.018 


(0.002) 


0.39 


(o 


031 


0.06 


(0.01) 


30 


W 


6D 4 

WW . t: 


44 8 

11 . CJ 


0.278 


(0.005) 


-0.018 


(0.004) 


0.52 


(o 


WZj J 


-0.06 


(0.02) 


O X 


1 


6D 8 

WW . (J 


^1 7 

W X . 1 


-0.534 


(0.005) 


0.534 


(0.009) 


0.54 


(o 


031 

WO I 


-0.37 


(0.04) 


;j X 


2 


61 3 


128 

X Zj W . VJ 


-0.534 


(0.005) 


0.534 


(0.009) 


0.49 


(o 


03") 

WO i 


-0.29 


(0.03) 


31 

O X 


W 


61 7 


81 2 

(J X . Zj 


-0.534 


(0.020) 


0.534 


(0.032) 


0.45 


(o 


WW 1 


-0.20 


(0.07) 


32 

OZj 


i 

X 


62 1 

WZj . X 


1 2Q 7 

X Zj t_/ . 1 


0.000 


(0.001) 


0.000 


(0.001) 


0.01 


(o 


01) 

VI) 


0.01 


(0.01) 


32 

OZj 


9 

Zj 


62 ^ 

WZj . W 


622 ^ 

WZjZj . W 


0.000 


(0.001) 


0.000 


(0.001) 


0.01 


(o 


011 

VI) 


-0.01 


(0.01) 


32 

OZj 


W 


62 Q 

WZj . <j 


Q18 1 

<J X(J . X 


0.000 


(0.000) 


0.000 


(0.000) 


0.00 


(o 


011 

VI) 


0.00 


(0.01) 


32 

OZj 


4 


63 4 

WW . t: 


^20 8 

W Zj W . O 


0.000 


(0.001) 


0.000 


(0.001) 


-0.01 


(o 


01) 

VI) 


0.01 


(0.01) 


33 


1 

X 


62 ^ 

WZj . W 


4D 7 

TlW . I 


-0.905 


(0.002) 


-1.579 


(0.006) 


0.34 


(o 


01 1 

VI) 


-0.16 


(0.01) 


33 


2 

Zj 


62 Q 

WZj . u 


60 Q 

WW . <J 


-0.905 


(0.004) 


-1.579 


(0.010) 


0.38 


(o 


011 

VI) 


-0.23 


(0.02) 


33 


O 


63 4 

WW . t: 


71 Q 


-0.905 


(0.009) 


-1.579 


(0.013) 


0.40 


(o 


02) 

WZj J 


-0.27 


(0.03) 


33 

■J -J 


4 


63 8 


52 7 

WZj . I 


-0.905 


(0.014) 


-1.579 


(0.021) 


0.38 


(o 


041 


-0.28 


(0.06) 


34 


1 


-28.2 


3.6 


-2.054 


(0.198) 


-8.488 


(0.243) 


-0.17 


(o 


41) 


-0.25 


(0.50) 


34 


2 


-27.8 


5.5 


-2.054 


(0.163) 


-8.488 


(0.208) 


-0.04 


(o 


34) 


-0.07 


(0.43) 


34 


3 


-25.7 


78.5 


-2.055 


(0.193) 


-8.489 


(0.182) 


-0.26 


(o 


42) 


-0.43 


(0.38) 


35 


1 


39.8 


6.0 


-2.766 


(0.026) 


-7.023 


(0.023) 


0.11 


(o 


09) 


-2.81 


(0.10) 


35 


2 


40.2 


15.3 


-2.767 


(0.022) 


-7.023 


(0.022) 


-0.06 


(o 


07) 


-2.96 


(0.05) 


35 


3 


40.6 


9.9 


-2.767 


(0.015) 


-7.023 


(0.012) 


-0.22 


(o 


04) 


-3.09 


(0.04) 


36 


1 


42.3 


0.7 


-2.765 


(0.017) 


-7.022 


(0.015) 


-0.07 


(o 


05) 


-2.91 


(0.05) 


36 


2 


42.7 


1.7 


-2.765 


(0.027) 


-7.022 


(0.012) 


-0.22 


(o 


09) 


-2.86 


(0.04) 


37 


1 


111.8 


107.2 


-2.787 


(0.034) 


-6.876 


(0.049) 


-0.70 


(o 


10) 


1.61 


(0.16) 



31 



Table 6 — Continued 



ID 


# 






x 




Vo 










(km s _1 ) 


(Jy/b) 


(") 




(") 


(mas yr -1 ) 


(mas yr -1 ) 


37 


2 


112.2 


59.5 


-2.787 (0.045) 


-6.; 


376 (0.051) 


-0.71 (0.11) 


1.56 (0.12) 


37 


3 


112.6 


16.8 


-2.787 (0.041) 


-6.; 


376 (0.049) 


-0.76 (0.09) 


1.57 (0.10) 


37 


4 


113.0 


6.1 


-2.787 (0.037) 


-6.; 


376 (0.037) 


-0.88 (0.08) 


1.67 (0.08) 



Note. — Individual H 2 maser spots in W51 Main/South that are grouped as features 
in Table [3j Column 1 indicates maser feature IDs corresponding to Table |3j Column 2 
shows maser spot IDs in different velocity channels in each feature. 



